The Cosmic Gravitational- Wave Background in a Cyclic Universe 
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Inflation predicts a primordial gravitational wave spectrum that is slightly "red," i.e. nearly 
scale-invariant with slowly increasing power at longer wavelengths. In this paper, we compute both 
the amplitude and spectral form of the primordial tensor spectrum predicted by cyclic/ekpyrotic 
models. The spectrum is exponentially suppressed compared to inflation on long wavelengths, and 
the strongest constraint emerges from the requirement that the energy density in gravitational waves 
should not exceed around 10 per cent of the energy density at the time of nucleosynthesis. 



The recently-proposed cyclic modelQ, 0] differs radi- 
cally from standard inflationary cosmology 0, Q , while 
retaining the inflationary predictions of homogeneity, 
flatness, and nearly scale-invariant density perturbations. 
It has been suggested that the cosmic gravitational wave 
background provides the best experimental means for dis- 
tinguishing the two models. Inflation predicts a nearly 
scale-invariant (slightly red) spectrum of primordial ten- 
sor perturbations, whereas the cyclic model predicts a 
blue spectrum. The difference arises because inflation 
involves an early phase of hyper-rapid cosmic accelera- 
tion, whereas the cyclic model does not. 

In this paper, we compute the gravitational wave spec- 
trum for cyclic models to obtain both the normalization 
and spectral shape as a function of model parameters, im- 
proving upon earlier heuristic estimates. We find that the 
spectrum is strongly blue. The amplitude is too small to 
be observed by currently proposed detectors on all scales. 
Hence, the discovery of a stochastic background of grav- 
itational waves would be evidence in favor of inflation, 
and would rule out the cyclic model. 

Readers unfamiliar with the cyclic model may consult 
[j| for an informal tour, and |(| for a recent analysis of 
phenomenological constraints. Cyclic cosmology draws 
strongly on earlier ideas associated with the "ekpyrotic 
universe" scenario. 0, 0, @ Briefly, the scenario can be 
described in terms of the periodic collision of orbifold 
planes moving in an extra spatial dimension, or, equiv- 
alently, in terms of a four-dimensional theory with an 
evolving (modulus) field cj> rolling back and forth in an 
effective potential V(d>). The field theory description is 
the long wavelength approximation to the brane picture 
in which the potential represents the interbrane interac- 
tion and the modulus field determines the distance be- 
tween branes. For the purposes of this paper, the field 
theoretic description is more useful. 

The potential (Fig. 1) is small and positive for large 4>, 
falling steeply negative at intermediate cj>, and increasing 
again for negative tj>. Each cycle consists of the following 
stages: (1) <j> large and decreasing: the universe expands 
at an accelerated rate as V(<f>) > acts as dark energy; 




end 



FIG. 1: Schematic of cyclic potential with numbers represent- 
ing the stages described in the text. To the left of (f> en d, where 
the scalar kinetic energy dominates, we approximate V with 
a Heaviside function, jumping to zero as shown by the dashed 
line. 



(2) <j> intermediate and decreasing: the universe is dom- 
inated by a combination of scalar kinetic and potential 
energy, leading to slow contraction and to the generation 
of fluctuations; (3) <j> negative and decreasing (beginning 
at conformal time T en d < 0): the generation of fluctua- 
tions ends, 4> rolls past 4> en d and, in the four-dimensional 
description, the universe contracts rapidly, dominated by 
scalar field kinetic energy, to the bounce (r = 0) at which 
matter and radiation are generated; (4) <f> increasing from 
minus infinity: the universe remains dominated by scalar 
field kinetic energy, which decreases rapidly compared to 
the radiation energy; (5) (f> large and increasing (begin- 
ning at T r > 0): the scalar field kinetic energy red-shifts 
to a negligible value and the universe begins the radiation 
dominated expanding phase; (6) <fi large and nearly sta- 
tionary: the universe undergoes the transitions to matter 
and dark energy domination, and the cycle begins anew. 
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We model the scalar field potential as: 

V(<f>) =V (1- e^/ M ?')<d(<f> - <t> end ) (1) 

where M v i is the reduced Planck mass and Q((j>) the 
Heaviside step function. A potential of this form, with an 
exponentially steep form, is required by the cyclic model 
in order to produce an acceptable spectrum of cosmolog- 
ical perturbations. 0, Q Choosing c = 10 for example 
results in a scalar spectral index n s — .96 which is com- 
patible with current constraints. The Heaviside function 
marks the end of the steeply decreasing part of the poten- 
tial; for <j) < <j) en d the potential is small and the universe 
is dominated by scalar field kinetic energy. 

Our calculation begins in the "ekpyrotic phase," stage 
(2), with the Einstein-frame scale factor contracting: 
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where a = 2/(c 2 — 2) -C 1 and r e fc = (1 — 2a)r en( j, being 
the conformal time the potential would have diverged to 
minus infinity had the exponential form continued. At 
t = T en d, the ekpyrotic phase ends and the "contracting 
kinetic phase," stage (3), begins: 
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T end < T < . (3) 



At t = 0, the universe bounces and the "expanding ki- 
netic phase," stage (4), begins: 
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Radiation is produced at the bounce, but is less than 
the scalar kinetic energy until, at r = r r , the expanding 
kinetic phase ends, and standard radiation-dominated, 
matter-dominated, and dark-energy-dominated epochs 
ensue. The transition times, r r and T en d, are given by 
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(5) 



(6) 



where H r = H(r r ) is the Hubble constant at r r , V en d — 
— V((f> en d) is the depth of the potential at its minimum, 
and % 1 is a small positive constant that measures the 
amount of radiation created at the bounce. Note that 
a(r) and a'(r) are both continuous at the transition time 
t = T en d, and we have chosen to normalize a(r) to unity 
at the start of radiation domination (a(r r ) = 1). 



THE PRIMORDIAL SPECTRUM, Ah(k,r r ) 

A quasi-stationary, isotropic, stochastic background 
of gravitational waves is characterized by the quantity 
Ah(k,r), the rms dimensionless strain per unit logarith- 
mic wavenumber at time r (i.e. the SL/L that would be 
measured by a detector with sensitivity band centered on 
mode k and bandwidth Afc = k). Accounting for both 
polarizations, it is given by Ah(k,r) — fc 3 / 2 |ft.fc(Y)|/7r, 
where the Fourier amplitude /ife(r) satisfies 



til + 2-h' k + k 2 h k = 
a " 



(7) 



and the boundary condition that the solution approaches 
the Minkowski vacuum at short distances 



h k (r) 
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To solve equation (JJJ, it is useful to define /fe(r) 
f( r ) hk{r) and rewrite as 



(8) 
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During the ekpyrotic phase, a(r) is given by |2|l. and the 
general solution of Q is 

f k (r) = ^(A 1 (k)H n 1 \y)+A 2 {k)H n 2 Hy)) , (10) 

where Ai^{k) are arbitrary constants, n = h — a, y = 

— k(r — Tek), an d H n ' are the Hankcl functions. The 
boundary condition (JSJ) implies 

Mk)= l - X [?r, A 2 (k)=0, 



(11) 



where we have dropped a physically irrelevant phase. In 
the contracting kinetic phase, stage (4), a(r) is given by 
(PJ), and the general solution of is 

f k (T) = V^(B 1 (k)H^ 1 \~kT) + B 2 (k)H^\-kT)j (12) 

where B\^{k) are arbitrary constants. Then, continuity 
of hk and h' k at r = r en d implies 
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(13) 



where x e = k\T en d\. Finally, in the expanding kinetic 
phase, a(r) is given by and the general solution of 
© is 

A(t) =Vkr(c 1 (k)H^(kr)+C 2 (k)H^(kr)) . (14) 
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To fix C\ 2(k), we need to match the solution across r = 
0. At the level of quantum field theory in curved space- 
time, the choice is essentially unique 11], and amounts 
to analytically continuing the positive (negative) fre- 
quency part of hk = fk/a around the origin in the lower 
(upper) half of the complex r-plane, so Hq 1 ' 2 ^— kr) — > 

r(2,l 



Hp '(kr). This yields 



(15) 



The pre-factor arises because a(r) differs by a factor of 
y/1 + x between the kinetic contracting and expanding 
phases; see Eqs. © and (@J. Combining our results, we 
arrive at the "primordial" dimensionless strain spectrum 
at the beginning of the radiation dominated epoch: 
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Ah(k,r r ) = (k 2 /nM pl ) v/2(l + X )i 



B 2 {k)H^\x r ) 
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(16) 



where x r = kr r and k < k en d- For k > k en d, the 
spectrum is cut off because these modes are not ampli- 
fied and, instead, Eq. i|16|) converges to the result for a 
static Minkowski background. Recall, in the string the- 
ory context, we are describing the collision of two orbifold 
planes and in the vicinity of the collision the background 
space-time becomes a nearly-trivial, flat, compactified 
Milne geometry, which is locally Minkowski. [8j After the 
bounce, in the Einstein-frame description we use here, 
the universe is in an expanding phase where modes no 
longer exit the horizon. Super-horizon modes are frozen, 
and sub-horizon modes are well described by the adia- 
batic (WKB) approximation so that new gravitational 
waves are not generated. 



THE PRESENT-DAY SPECTRUM, Ah(k,T ) 

To convert from the primordial spectrum to the 
present-day spectrum Ah(k,Ta) = Th(k) Ah(k , r r ) we 
need to know the transfer function, Th(k). To ap- 
proximate Th{k), note that Ah(k,r) is roughly time- 
independent outside the horizon, and decays as a -1 once 
a mode re-enters the horizon. Therefore, the transfer 
function is ~ 1/(1 + z r ) for modes already inside the 
horizon at the onset of radiation domination r r , and 
~ 1/(1 + Zk) for modes that entered at red shift Zk be- 
tween r r and To. Using the fact that H cx a~ 2 dur- 
ing radiation domination and H cx a~ 3 / 2 during matter 
domination (neglecting the change in <?*) we find 
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where fco = o-qHq, k eq = a eq H eq , k r = a r H ri and 
k enc i = a en d\H en d\ denote the modes on the horizon to- 
day (to), at matter-radiation equality (r eg ), at the start 



FIG. 2: A schematic comparison of the dimensionless strain 
observed today Ah(k,ro), as predicted by inflation and the 
cyclic model. Here tit is the inflationary tensor spectral index 
(a small negative number), and a < 1 in the cyclic model is 
a small positive number. k r denotes the mode on the horizon 
at the start of radiation domination. 



of radiation domination (r r ), and at the end of the ekpy- 
rotic phase (r en d), respectively. 

The gravitational wave spectrum can be divided into 
three regimes. There is a low frequency (LF) regime cor- 
responding to long wavelength modes that re-enter after 
matter-radiation equality (k < k eq ), and a medium fre- 
quency (MF) regime consisting of modes which re-enter 
between equality and the onset of radiation domination 
(k eq < k < k r ). (We ignore the recent dark energy dom- 
inated phase, which has negligible effect.) The spectrum 
for these two regimes is: 



Ah. 



r*fcg 



k -l+a 

k a /k eq 



(LF) 
(MF) 



(18) 



Finally, modes which exit the horizon during the ekpy- 
rotic phase (before T en d), and re-enter during the expand- 
ing kinetic phase (after the bound but before r r ) result 
in a high frequency (HF) band (k r < k < k en d)' 
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The HF band runs over a range k en d/k r = T, and this 
quantity is strongly constrained by the requirement that 
the scalar field cross the negative region of the poten- 
tial before radiation domination begins, which requires 
that 6] 
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where V$ is today's value of the dark energy density. This 
equation, combined with ©, gives a lower bound on T, 

. For example, for V en d around the 
GUT scale and c = 10, we find T > 10 8 . Fig. H schemat- 
ically depicts Ah(k,To) in the cyclic scenario and com- 
pares it to the inflationary spectrum, [l^ 

Another useful quantity is Q gw (k, tq), the gravitational 
wave energy per unit logarithmic wavenumber, in units 
of the critical density |lfj, [l2( : 



7-*W = U& A " lW " ,2:l! 
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In the cyclic model, fl gw (k 1 To) is very blue, with nearly 
all the gravitational wave energy concentrated at the 
high-frequency end of the distribution. 
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OBSERVATIONAL CONSTRAINTS AND 
DETECTABILITY 

The strongest observational constraint on the gravi- 
tational spectrum in the cyclic model comes from the 
requirement that the successful predictions of big bang 
nucleosynthesis (BBN) not be affected, which requires 
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From the above equations, (|19ll and l|21|l . and using 

i i 

1 + z eq s» k1 q /k,Q, and T r ~ Hr Mp t for the temperature 
at radiation domination, we obtain a total £1 in gravi- 
tational waves of ~ (2aV end /T r M Pl )^ [36tt 3 (1 + z e9 )]~\ 
which from (|22|l implies 



(23) 



where, for simplicity, we have ignored the factor which 
depends on the number of thermal degrees of freedom, 
which further weaken this bound. 

The other observational constraints are much weaker. 
[l3| From the CMB anisotropy, one infers Ah{f ~ 
10 _18 Hz) < 10 -5 ; from precision pulsar timing, Ah(f ~ 
10- 8 Hz) < 10" 14 . Optimistic goals for LISA and 
advanced LIGO are strain sensitivities of Ah(f ~ 
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spectively. Fig. 3 shows results for values of T r and V en d 
consistent with all constraints on the cyclic model 0. 

Even if the parameters are chosen to saturate the BBN 
constraint, the spectrum is still orders of magnitude be- 
low the sensitivity of anticipated instruments. In partic- 
ular, it is hopeless to search for an imprint in the polar- 
ization of the CMB from cyclic model gravity waves be- 
cause the predicted amplitude on large scales is so small. 
Hence, the detection of a stochastic gravitational wave 
imprint in the CMB polarization would be consistent 
with inflation and definitively rule out the cyclic model. 



FIG. 3: The present-day dimensionless strain, Ah(k,ro), pre- 
dicted by the cyclic model with T r = 10 7 GeV and = 
10 14 GeV. These parameters yield a gravity wave density four 
orders of magnitude below the BBN bound. Some observa- 
tional bounds and (optimistic) future strain sensitivities are 
indicated. The dashed arrows mean the empirical bounds lie 
well above range of Aft displayed here. 



We thank A. Tolley and J. Khoury for helpful conversa- 
tions. LB is supported by an NSF Graduate Fellowship. 
This work was also supported in part by US Depart- 
ment of Energy grant DE-FG02-91ER40671 (PJS) and 
by PPARC-UK (NT). 



[1] 

[2] 

[3 
[4 

[5] 
[6] 

[7 

[8 

[9 

[10 

[11 



[12 
[13 



P. J. Steinhardt and N. Turok, Science 296, 1436-1439 
(2002). 

P. J. Steinhardt and N. Turok, Phys. Rev. D 65, 126003 
(2002). 

A. H. Guth, Phys. Rev. D 23, 347 (1981). 
A. D. Linde, Phys. Lett. B 108, 389 (1982); A. Albrecht 
and P. J. Steinhardt, Phys. Re v. Lett. 48, 1220 (1 982). 
P. J. Steinhardt and N. Turok, astro-ph/0204479 

J. Khoury^ P J. Steinhardt, and rT! ^Turok, 

hep-th/0307132 2003. 

j. Khoury, B. A. Ovrut, P. J. Steinhardt, and N. Turok, 
Phys. Rev. D 65, 123522 (2002). 

J. Khoury, B. A. Ovrut, N. Seiberg, P. J. Steinhardt, and 
N. Turok, Phys. Rev. D 65, 086007 (2002). 
J. Khoury, B. A. Ovrut, P. J. Steinhardt, and N. Turok, 
Phys. Rev. D 66, 046005 (2002). 

K. S. Thorne, in 300 Years of Gravitation, edited by S. 
Hawking and W. Israel (Cambridge: Cambridge Univer- 
sity Press, Cambridge, England, 1987), pp. 330-458. 

A. J. Tolley and N. Turok, Phys. Rev. D 66, 106005 
(2002); A.J. T olley, N. Turok and P.J. Steinhardt, 
hep-th/0306109 (2003). 

M. S. 'lYi rneiTPhys. Rev . D 55, 435 (1997). 

B. Allen, gr-qc/9604033 K. S. Thorne, gr-qc/9506086 



